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SUMMARY 

This  grant  provided  resources  to  study  (1)  ionospheric  disturbances  that  cause  severe  disruptions  of 
radio  communications  and  navigational  systems  crucial  to  Department  of  Defense  programs,  and  (2) 
tsunami-generated  ionospheric  signatures  of  waves  that  can  be  used  to  study  both  upward  and 
horizontal  ocean -atmosphere  coupling.  Our  primary  observational  technique  is  optical — one  that 
provides  broad  coverage  of  upper  atmosphere  structures  and  disturbances.  Using  instrumentation 
provided  by  the  Defense  University  Research  Instrumentation  Program  (DURIP)  we  conducted  our 
studies  from  several  sites  in  North  and  South  America.  Each  site  was  equipped  with  a  Boston 
University  state-of-the-art  optical  instrument  called  an  AU-Sky-magier  (ASI),  meaning  that  it  records 
atmospheric  emissions  using  a  180°  field-of-view  achieved  by  a  “fish-eye”  lens.  Moreover,  we  located 
these  instruments  at  sites  that  share  a  common  geomagnetic  field  line — called  conjugate  points.  For 
example,  the  ASI  at  the  Arecibo  Radio  Observatory  in  Puerto  Rico  has  its  geomagnetic  conjugate  point 
in  Mercedes  (Argentina),  while  our  ASI  at  the  higher  latitude  Millstone  Uill/Haystack  Observatory 
(MA)  has  its  conjugate  point  covered  by  an  ASI  in  Rothera  (A\ntartica). 

I .  Ionospheric  disturbances  at  low  latitudes.Of  particular  interest  to  DoD  radio  communications 
and  navigational  systems  are  the  radio  disruptions  ( signal  amplitude  and  phase  scintillations )  produced 
by  ionospheric  irregularities  at  equatorial  and  low  latitudes  from  sunset  to  dawn.  These  are  called 
Equatorial  Spread-F  (ESF)  scintillations  and  they  are  the  strongest  anywhere  on  the  globe.  The  optical 
signatures  of  such  irregularities  can  be  captured  by  ASI  methods  that  span  thousands  of  kilometers. 

We  use  the  oxygen  ‘red  line’  at  6300  A  to  make  images  of  where  ionospheric  “airglow  depletions” 
exist  at  times  throughout  the  night.  Any  radio  propagation  ray  path  (e.g.,  using  GPS)  through  one  of 
these  depletions  can  experience  strong  scintillation  to  the  point  of  total  signal  disruption  and  therefore 
the  loss  of  an  important  navigational  tool.  Such  overall  scenarios  are  called  Space  Weather  events. 

Starting  in  October  2014,  the  BU  ASI  network  established  the  first-ever  set  of  three  all-sky-  imagers 
spanning  latitudes  from  the  geomagnetic  equator  to  low  latitudes  in  each  hemisphere.  The  three  ASI 
observatories  are  located  in  three  South  American  countries  linked  by  a  common  set  of  geomagnetic 
field  lines:  Villa  de  Leyva  (Colombia)-Jicarmarca  (Peru)-EI  Leoncito  (Argentina).  In  Figures  1  and  2, 
we  show'  how  this  new  resource  can  be  used  on  a  given  night  (October  30,  2014)  to  study  the  latitudc- 
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altitude  relationship  between  6300  A  airglow  depletion  signatures  and  the  ionospheric  irregularities 
that  cause  radio  disruptions.  The  key  features  to  note  are  described  in  the  captions.  We  stress  that 
using  single-site  images  from  any  one  of  these  three  stations  it  would  be  impossible  to  describe 
accurately  the  spatial-temporal  evolution  of  ESF  disturbances  on  this  night.  Moreover,  the  coherence 
of  the  perturbations  effects  can  be  used  to  specify  space  weather  environmental  impacts  in  area-denied 
locations  (e.g»  information  from  Argentina  can  tell  if  the  ionosphere  above  Columbia  is  experiencing 
radio  disruptions).  This  type  of  research,  enabled  by  this  ONR  grant,  will  be  conducted  by  new  grant 
funds. 
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Figure  1.  First  example  of  3-site,  all -sky- imager 
depictions  of  6300  A  airglow  depletions 
associated  with  trans-equalorial  plasma 
instabilities.  Two-minute  exposure  images  taken 
“simultaneously”  during  the  pre-midnight  hours 
of  30  October  2014  were  achieved  using 
instruments  at  Villa  dc  Leyva  (Colombia), 
Jicamarca  (Peru),  and  FJ  Leone i to  (Argentina).  In 
the  central  image,  the  N-S  aligned  airglow 
depletion  through  zenith  shows  where  ionospheric 
irregularities  are  located  in  regions  close  to  the 
geomagnetic  equator.  This  dark  feature 
“connects”  to  structured  airglow  depletions  at 
higher  latitudes  to  the  north  and  south,  indicating 
that  an  entire  magnetic  meridian  experiences 
irregularities.  The  magnetic  conjugacy  is  evident 
with  depletions  extending  to  ±  20°  magnetic 
latitudes  north  and  south  of  the  magnetic  equator. 
A  similar  pair  of  conjugate  airglow  depletions  are 
captured  to  the  east,  but  its  equatorial  signature 
falls  beyond  the  fleld-of-view  of  the  equatorial 
station.  The  degree  of  structuring  appears  more 
pronounced  in  the  northern  hemisphere. 
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Figure  2.  Using  the  same  format  as  in  Figure  I , 
these  three  images  are  for  the  post-midnight 
period  on  the  same  night  of  30  October  20 1 4.  A 
dark  airglow  depletion  spans  the  image  from  the 
equatorial  site  (Jicarmarca).  Ionospheric 
irregularities  “visualized”  by  their  airglow 
depletions  are  far  less  evident  at  sites  north  and 
south  of  the  geomagnetic  equator.  There  is  a  hint 
of  highly-tilted  depletions  at  the  western  edge  of 
the  FOVs  at  Villa  de  Leyva,  not  seen  at  El 
Leoncito,  The  complex  airglow  depletions  to  the 
east  at  El  Leoncito  do  not  have  an  easily- 
identified  presenee  in  the  images  from  Villa  de 
Leyva.  Events  such  as  depicted  in  Figures  I  and  2 
are  now  capable  of  being  analyzed  systematically 
for  the  first  time. 


2.  Studies  of  Earthquake  and  Tsunami-induced  Waves  in  the  Ionosphere . 

One  of  the  more  spectacular  uses  of  all-sky  airglow  imaging  was  the  recent  discovery  of  waves  in  the 
ionospheric  airglow  layer  caused  by  the  great  earthquake  and  tsunami  of  1 1  March  201 1  (Makela  el  ah, 
201 1).  They  described  a  large  thermospheric  gravity  wave  captured  in  6300  A  images  that  propagated 
over  the  Hawaiian  archipelago  almost  five  hours  after  the  occurrence  of  the  earthquake. 

Most  interesting,  Makela  et  al.  (201 1)  reported  that  the  airglow  wave  signatures  over  Hawaii  arrived 
approximately  an  hour  anead  of  the  tsunami — thus  suggestive  of  a  possible  forecast  or  early-warning 
capability  from  AS  I  observations. 

As  shown  in  Figure  1 ,  Boston  University  operates  an  ASI  at  the  El  Leoncito  Observatory  in  Argentina 
and  Figure  3  displays  three  red-line  images  from  El  Leoncito  showing  the  passage  of  gravity  wave 
signatures  on  the  night  of  12  March  201  L  The  analyses  of  such  images  yielded  a  horizontal  phase 
speed  of  204  ±  4  m/sec,  a  wavelength  of  200  ±  23  km.  and  period  of  1 6  ±  2  minutes. 

The  great  circle  ground-distance  between  the  El  Leoncito  Observatory  and  the  earthquake  epicenter 
is!  7,750  km  for  airglow  layer  at  250  km  altitude  (Figure  4).  Using  the  measured  arrival  time  at  El 
Leoncito  and  horizontal  phase  speed  of  the  waves,  the  estimated  “start-time”  of  the  tsunami-genic 
gravity  wraves  in  Japan  came  out  to  be  05:46:24  UT  (±5,5  minutes)  on  1 1  March.  This  is  within  1.5 
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minutes  of  the  recorded  onset  time  in  Japan,  well  within  the  uncertainties  associated  with  image 
processing.  The  waves  detected  at  El  Leone ito  thus  travelled  for  24 hr  8min  ±5  min  to  nearly  the 
opposite  side  of  the  globe. 


Figure  5  shows  tsunami  arrival  time  data  from  a  series  of  tidal  monitoring  stations  located  along  the 
Chilean  coast,  as  recorded  by  the  NOAA's  National  Geophysical  Data  Center 

[http://www.ngdc.noaa.gov/hazard/].  It  is  clear  from  this  figure  that  there  is  no  significant  difference  in 
the  arrival  times  of  the  airglow  and  ocean  signatures  of  the  disturbance.  This  is  in  marked  contrast  to 
the  Makela  et  al.  (201 1 )  results  showing  the  airglow  waves  over  Hawaii  arriving  about  one  hour  ahead 
of  the  tsunami  waves.  The  reason  for  such  differences  are  unclear,  and  thus  a  topic  in  need  of  further 
observational  and  modeling  efforts.  One  possibility  could  be  that  variations  in  the  propagation  speed  of 
the  tsunami  result  from  variations  in  ocean  depth  and/or  temperature  between  vastly  different  regions 
of  the  Pacific  Ocean  (Smith  et  al.,  2015). 


Figure  3.  Sequence  of  three  TD  all- 
sky  images  at  the  El  Leone  ito 
Observatory  on  12  March  205  1.  The 
unwarped  630.0  nm  images  show  the 
south-eastward  progression  ol'  an 
extensive  thermospheric  gravity  wave 
event.  The  schematic  tracings  in  the 
lower  part  of  each  panel  highlight  the 
waves  more  clearly. 


Figure  4.  Global  map  showing  the 
relative  locations  of  the  earthquake 
epicenter  and  the  El  Leoncito 
Observatory  (ELO)  in  Argentina.  The 
propagation  path  of  the  tsunami  to  El 
Leoncito,  a  distance  of  almost  1 8.000 
km,  is  shown.  The  propagation  path  to 
Hawaii  (viz-a-viz  Makela  et  al. 

(20 1  I ))  is  also  shown  for  comparison. 
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Figure  5.  The  measured  oceanic  tsunami 
travel  times  to  the  NOAA  Chilean 
coastal  tidal  monitoring  stations  (data 
from  NDGC)  with  their  great  circle 
distances  from  the  earthquake  epicenter. 
The  zenith  crossing  time  of  the  leading 
630.0  nm  wavefront  at  El  Leone i to 
(ELO)  is  also  plotted. 


Program  matics 

The  success  of  this  grant  was  made  possible  by  an  existing  research  infrastructure  within  the  Imaging 
Science  Laboratory  of  the  Center  for  Space  Physics  at  Boston  University. 

Key  Individuals 


(a)  Michael  Mendillo,  Professor  of  Astronomy,  served  as  PI  and  directed  the  overall  programs 
of  observations  and  data  analysis. 

(b)  Jeffrey  Baumgardner,  Senior  Research  Scientist  in  the  Center  for  Space  Physics,  designs, 
constructs  and  repairs  all  instrumentation;  he  participates  in  data  analysis  and  interpretation. 

(c)  Joei  Wroten,  Senior  Staff  Researcher,  is  in  charge  of  data  analysis  and  archiving;  she 
maintains  our  website,  conducts  image  processing,  and  works  with  the  PI  on  ionospheric 
effects  associated  with  solar  activity. 

(d)  Carlos  Martinis  (Research  Assistant  Professor)  leads  the  analysis  and  interpretation  of  the 
imaging  data  from  optical  imaging  systems  at  equatorial  and  low  latitudes. 

(e)  Dr.  Steven  Smith  (Senior  Research  Scientist)  leads  the  analysis  and  interpretation  of  wave 
disturbances  in  the  upper  atmosphere. 

(f)  During  the  three-year  period  covered  by  this  grant,  several  undergraduate  research  assistants 
worked  with  the  senior  staff  members  listed  above — leaning  data  analysis  methods  and 
developing  skills  in  science  and  engineering  linked  to  space  research. 
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WORK  COMPLETER 


We  have  shown  that  multi-station  observations  at  conjugate  points  of  the  geomagnetic  field  can  be 
used  to  characterize  the  spatial  and  temporal  patterns  of  ionospheric  irregularities  that  cause 
communications  disruptions.  In  addition,  we  have  demonstrated  that  tsunami  induced  waves  in  the 
upper  atmosphere  are  capable  of  propagation  to  distances  well  beyond  previous  estimates.  The 
question  of  using  such  waves  to  predict  tsunami  arrival  times  is  one  in  need  of  additional  research. 

IMPACT/APPLICATIONS 

We  have  conducted  a  pilot  study  that  relates  all-sky-imaging  in  one  hemisphere  to  that  in  another — 
with  the  images  showing  where  ionospheric-disruption  effects  would  occur.  We  found  a  broad 
agreement  on  the  time  scale  of  a  night,  but  a  lack  of  coherence  in  space  and  time  on  time  scales  of  an 
hour  or  less.  This  implies  that  plasma  instabilities  may  be  affected  by  local  conditions  more  than 
previously  considered. 

RELATED  PROJECTS:  Synergistic  Use  of  Data 

At  several  o:  the  observatories  where  BU  all-sky-imagers  operate  to  gather  data  about  the  ionosphere, 
our  background  (“control  wavelength”)  data  reveal  where  clouds  are  located.  This  is  of  use  to 
colleagues  operating  telescopes  on  the  same  mountain,  or  by  remote  control  from  home  institutions,  as 
a  reliable  source  of  “seeing”  conditions  for  their  experiments  of  astronomical  targets. 

CONCLUSIONS 

The  use  of  our  DURIP  and  ONR  sponsored  research  instruments  to  support  the  highly-focus  research 
goals  described  in  this  report  represent  a  state-of-the-art  usage  of  DoD  research  funds.  They  have  a 
clear  applications  avenue  for  space  weather  now-casting  and  predictions.  Our  studies  also  contribute 
in  fundamental  ways  to  the  advancement  of  science  related  to  the  broad  field  geo-space  environment. 
For  the  training  of  new  scientist  and  engineers  enabled  by  this  grants,  wc  express  our  gratitude  to  the 
Office  of  Naval  Research  and  its  program  managers.  From  an  international  relations  perspective,  this 
ONR/BU  program  also  had  a  very  positive  set  of  interactions  with  scientific  and  educational 
institutions  in  South  America. 

Finally,  there  were  no  inventions  or  patents  produced  from  the  research  supported  by  this  grant. 

REFERENCES 

Makela,  J.  el  al.  (201 1),  Imaging  and  modeling  the  ionospheric  airglow  response  over  1  lawaii  to  the 
tsunami  generated  by  the  Tohoku  earthquake  of  1 1  March  2011,  Geophys.  Res.  Lett.,  L20807, 
doi:l  0.1 029/20 1 IGL04760. 

Smith,  S.,  C.  Martinis,  J.  Baumgardner,  and  M.  Mendillo,  All-sky  imaging  of  transgiobal 
thermospheric  gravity  waves  generated  by  the  March  1 1  Tohoku  Earthquake  (2015),  .1.  Geophys.  Res. 
Space  Physics,  120,  doi:t 0.1002/20 15JA02 1 638. 


6 


PUBLICATIONS  and  PRESENTATIONS  RESULTING  PROM  THIS  GRANT  (all  or  in  part) 
—  entries  listed  year — BU  authors  in  bold  type — 

2014 

Pallamraju,  D..  J.  Baumgardner  &  8  co-authors,  Daytime  wave  characteristics  in  the 
mesosphere  and  lower  thermosphere  region:  Results  of  the  Ballon-borne 
Investigations  of  Regional-atmospheric  Dynamics  Experiments,  J.  Geophys.  Res., 
doi :  1 0 . 1 002/20 1  3.1  AO  1968,  2014. 

Makela.  J.,  and  12  co-authors  including  C.  Martinis,  Storm  time  response  of  the 
midlatitude  thermosphere:  Observations  from  a  network  of  Fabry-Perot 
interferometers,  .1.  Geophys.  Res.,  1  1 9.  8  doi:  1 0. 1 002/20 1 4ja0 1 9832,  2014. 

Martinis,  C.  et  al.  (including  R.  Macinnis,  C.  Sullivan,  J.  Baumgardner,  and  M. 

Mendillo),  Ionospheric  studies  from  low  to  sub-auroral  latitudes  using 
geomagnetic  conjugate  all-sky  imagers  and  satellite  data.  Airier.  Geophys.  Union 
meeting,  San  Francisco,  December  2014. 

Harding,  et  al.  (including  C.  Martinis),  Anomalous  Doppler  Shift  in  the  storm-time 
middle  red-line  emission,  Amer.  Geophys.  Union,  San  Francisco,  Dec.,  2014, 

Gentile  et  al.  (including  C.  Martinis),  C/NOFS  and  the  Air  Force  Research  Laboratory’s 
conjugate  point  study,  Amer.  Geophys.  Union,  San  Francisco,  Dec.,  2014. 

Smith,  S.M.,  Optical  Mesospheric  Instruments  During  Deepwave.  Deepwave  Planning  Meeting, 
Boulder,  CO.,  2014. 

Fritts,  D.,  R.  Smith,  M.  Taylor,  .1.  Doyle,  S.  Eckcnnann,  and  S.  Smith,  The  Deep  Propagating  Gravity 
Wave  Experiment  (DEEPWAVE):  Science  Overview  and  Approach,  Deepwave  Science  and 
Operations  Meeting,  Christchurch,  N.Z.,  21-23  Jan  2014. 


2015 

Stefanello,  M,,  M.  Muella,  D.  Amorim,  C  Machado.  J.  Bageston,  A.  Pimenta,  C.  Martinis,  and  N. 
Schuch,  01  630.0  nm  all-sky  image  observations  of  medium-scale  traveling  ionospheric 
disturbances  at  geomagnetic  conjugate  points.  J.  Atmos.  Solar-Terr.  Phys.,  128 .  58-69.  2015. 

Park,  J.  et  al.  (including  C.  Martinis),  A  dayside  plasma  depletion  observed  at  midlatitudes  during 
quiet  geomagnetic  conditions,  Geophys.  Res.  Lett.,  42(4),  967-974, 
doi:  1 0.1 002/20 14GL062655.  2015. 

Rodrigues,  F.,  et  al.  (including  C.  Martinis),  AMISR-14  observations  of  equatorial  spread  F,  Geophys. 
Res.  Lett..  42(  1 3),  5  1 00-5 1 08.  doi :  1 0. 1 002/20 1 5GL064574,  2015. 

Rodrigues,  F.,  M.  N  ico  I  Is,  M.  Milla,  J.  Smith,  R.  Varney,  A.  Strom  me,  J. 

Arratia,  and  C.  Martinis,  AMISR-14:  Observations  of  equatorial 


7 


spread  F,  Geophys.  Res.  Lett.,  42{  13),  5100-5108.  2015. 


Hickey,  D.,  C.  Martinis  et  al..  Concurrent  observations  at  the  magnetic  equator  of 
small-scale  irregularities  and  large-scale  depletions  associated  with 
equatorial  spread  F,  J.  Geophys.  Res.  Space  Phys.,  120  (12),  10,883- 
1 0,896,  doi :  1 0. 1 002/20 1 5JA02 1 99 1 ,  20 1 5 . 

Martinis,  C.,  J.  Baumgardner,  M.  Mcndillo,  J.  Wroten,  A.  Coster,  and  L.  Paxton,  The 
night  when  the  auroral  and  equatorial  ionospheres  converged,  J.  Geophys.  Res. 
Space  Physics,  1 20(9),  808508095,  doi:  1 0. 1 002/20 1 5J A02 1 555,  20 1 5. 


Martinis,  C.,  J.  Baumgardner,  M.  Mcndillo,  J.  Wroten,  R.  Mucin n is.  C.  Sullivan,  D. 
Hickey,  Conjugate  Behavior  of  MSTIDs  and  ESF  Structures,  AGU  Meeting,  San 
Francisco,  December,  2015. 

Wroten,  J.,  J.  Baumgardner,  M.  Mcndillo,  SAR  Arcs  We  Have  Seen:  Variability  of 
Ring  Current  -  Ionosphere  Interactions,  AGU  Meeting,  San  Francisco, 

December,  2015. 


2016 

Mendillo,  M.,  J.  Baumgardner,  and  J.  Wroten,  SAR  arcs  we  have  seen:  Evidence  for 
variability  in  stable  auroral  red  arcs,  J.  Geophys.  Res.  Space  Physics,  120, 
doi:  1_ 9.  W02/20UJA02 !_  72Z  2016. 

Burke,  W.  C.  Martinis,  P.  Lai,  L.  Gentile,  C,  Sullivan,  and  R.  Pfaff,  C/NOFS 

observations  of  electromagnetic  coupling  between  magnetically  conjugate 
MSTID  structures,  J.  Geophys.  Res.  Space  Phys.,  121(3),  2569-2582, 
doi :  1 0.1 002/20 15JA02 1965,  2016. 

Park,  J.,  C.  Martinis,  II.  Luhr,  R.  Pfaff,  and  Y.-S.  Kwak,  Hemispheric 
asymmetry  in  transition  from  equatorial  plasma  bubble  to  blob  as 
deduced  from  630.0nm  airglow  observations  al  low  latitudes,  J. 

Geophys.  Res.  Space  Phys.,  121(1),  881-893, 
doi:  1 0. 1 002/201 5J A022 1 75,  20 1 6. 

Burke,  W.,  R.  Pfaff,  C.  Martinis,  and  L.  Gentile,  C/NOFS  remote  sensing  of 
ionospheric  reflectance,  J.  Geophys.  Res.  Space  Phys.,  121(5), 

4924-4932,  doi :  10.1 002/20 16JA022345,  2016. 


Presentations  on  Synergistic  Topics:  Comparative  Ionospheres 

Mendillo,  M.,  Narvaez,  C.,  Lawler,  G.,  Kofman,  W.,  Mouginot,  .1.,  Morgan,  D.,  & 
Gurnett,  The  equivalent  slab  thickness  of  Mars’  ionosphere:  Implications 
for  thermospheric  temperature  Geophys.  Res.  Lett.,  42,  3560-3568, 
doi:  1 0. 1 002/20 1 5GL063096, 20 1 5. 


8 


